Ripening of Porous Media 



(N 
O 

o 

(N 

» 

o 
O 

(N 
(N 



43 

o 



CO 



I 

C 

O 

o 



> 
m 
O 
in 
o 

(N 
O 

c3 



13 

o 
o 



X 
S3 



Benny Davidovitch, Deniz Erta§ and Thomas C. Halsey 
Corporate Strategic Research, ExxonMobil Research and Engineering, 1545 Route 22 East, Annandale, NJ 08801 

(Dated: October 22, 2002) 

We address the surface tension-driven dynamics of porous media in nearly saturated pore-space 
solutions. We linearize this dynamics in the reaction-limited regime near its fixed points - surfaces of 
constant mean curvature (CMC surfaces). We prove that the only stable interface for this dynamics 
is the plane, and estimate the time scale for a CMC surface to become unstable. We also discuss 
the differences between open and closed system dynamics, pointing out the unlikelihood that CMC 
surfaces are ever realized in these systems on any time scale. 
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When solid grains are suspended in a solution satu- 
rated with the molecular constituents of the grains, they 
undergo coarsening under the thermodynamic driving 
force of surface tension. During this phenomenon, known 
as Ostwald ripening, the free energy of the system is low- 
ered by minimizing the contact area between the coexist- 
ing phases. Molecules dissolve from high-curvature areas 
of the interface, pass through the solution, and can pre- 
cipitate in low-curvature surface regions. This dynamics 
can be very complicated and depends on many parame- 
ters such as chemical composition, induced temperature 
and pressure fields, and other factors. 

In this Letter we study Ostwald ripening in porous me- 
dia, where the fluid in the pore space is approximately 
saturated with the ingredients of the solid phase. An 
example would be a sedimentary material such as sand- 
stone, with the water in the pore space saturated with the 
silica components of the rock. Another example would 
be crushed ice, with water vapor saturating the air in 
the pore space between ice grains. We identify the fixed 
points of dissolution-precipitation dynamics as surfaces 
of "Constant Mean Curvature", (CMC surfaces), and 
then show quite generally that these surfaces are unstable 
fixed points of this dynamics. 

If one suspends grains in a solution, a mean field the- 
ory due to Lifshitz-Slyozov and Wagner is successful in 
capturing the dynamics at late stages of the ripening pro- 
cess for low solid volume fraction Another example 
of surface tension driven ripening arises in the kinetics 
of foams |^| . Unlike these examples, in typical porous 
media both the solid and pore space components of the 
medium are connected. 

We thus pursue an approach more suited to porous 
media. Following || we consider the evolution of an 
interface T(t) = x(u(t),v(t)) between a porous, single 
component, isotropic solid and its ideal solution in the 
interstitial fluid. The solid is subject to a first order 
dissolution-precipitation reaction in a flow field of veloc- 
ity v. The normal velocity of the surface u n (x) into the 
pore space is given by 



where Kf is the dissolution rate, and the precipitation 
rate is controlled by the Boltzmann factor associated with 
the difference A/i(x) = /i sur (x) — /x so i(x) between the 
chemical potentials of solid and dissolved molecules at 
the interface T. Referred to the chemical potential /ifl at 
for a flat surface in equilibrium with an ideal saturated 
solution of concentration c sa t , these are given by 
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where v m is the molecular volume in the solid, er is the in- 
tcrfacial energy, c(x) is the concentration near the surface 
point x, and i?(x) = (1/2)(SS/5V) is the mean curva- 
ture, which measures the local variation in surface area 
SS with respect to a volume change SV of the solid. Here 
and elsewhere we define all concentrations with respect 
to the concentration in the solid. For A/z <C kT, Eq. (|l|) 
reduces to 
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Since H should be no larger than the inverse of a typical 
pore size L, H ~ L~ l , and v m a /kT is usually the inverse 
of a molecular scale, 
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we see that the linearized form of the dynamics holds for 
a nearly saturated solution c ~ c sa t ■ 

The transport of the concentration field c in the so- 
lution is described by the advection-diffusion equation, 
with the appropriate boundary condition at the solid in- 
terface ||. The velocity field satisfies the Navier-Stokes 
equation, with appropriate boundary conditions at infin- 
ity and at the solid surfaces inside the porous medium. 

Since the full solution of these equations in a porous 
medium is daunting, we concentrate on two limiting cases 
of this problem, which we term "perfect open" and "per- 
fect closed" systems. In a perfect open system, the ve- 
locity of the flow through the medium is considered suf- 
ficiently rapid that in effect each element of the surface 
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is in contact with fluid whose concentration c of the so- 
lute is fixed by a distant reservoir. This criterion can be 
achieved by taking the limit of small Damkohlcr num- 
ber Da = Kf/v — > 0, with v the average velocity of 
flow through the system 0. Most of this study will 
concentrate on this limit. The other limit of a "perfect 
closed" system, in which v — and the diffusion constant 
D — > oo, will be explored at the end of this Letter. 

In the perfect open case, the concentration c(x) on the 
solution side of the interface has the constant value Coo, 
and the interface dynamics is governed only by Eq. (Q), 
which yields 
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where = [kT / (2i/ m a)] log(coo / 'c sat ) . 

CMC surfaces, for which H(x) — H* everywhere on 
the surface, are the fixed points of the dynamics (||). 
These surfaces have been studied in a variety of contexts, 
notably for their relationship to certain phases of block 
co-polymers ||. The simplest examples are a sphere of 
radius r = jj and a cylinder of radius r = 

A special case of CMC surfaces are "minimal surfaces" , 
for which H(x) = 0. These surfaces have been exten- 
sively studied in the context of analytic function theory, 
and many such surfaces have been discovered Q. For 
= 0, Eq. (H) is the equation of motion for a mani- 
fold seeking to minimize its own surface dynamics 
that is sometimes called "motion under mean curvature" 
0. Of particular interest in the context of porous me- 
dia are the triply periodic minimal surfaces, of which the 
Schwartz P-surface, shown in Fig. 0(a), is a classic exam- 
ple. Note that both sides of the surface form connected 
components, as in most porous media. Anderson et al. 
have extended several such surfaces into CMC surfaces 
with H 7^ . The extension of the Schwartz P-surface 
to a family of CMC surfaces is shown in Fig. ||, with 
an example depicted in Fig. |](b). These surfaces vary 
continuously from the P-surface (M) to a simple cubic 
lattice of barely touching spheres at one endpoint (A), 
and to a simple cubic lattice of barely touching spherical 
holes at the other endpoint (A') (see Fig. ||). Anderson 
et al. also generated families corresponding to fee, bec, 
and diamond lattices of spheres, with qualitatively sim- 
ilar properties. Although we are aware of no reported 
examples of "amorphous" families of CMC surfaces, it is 
likely that such families exist [|| . Such families would be 
most relevant to real porous media. 

We now show that all CMC surfaces with the exception 
of planes are unstable fixed points of Eq.(||). Consider a 
surface x'(u, v) = x + e(x)n(x), whose deviation from a 
CMC surface x(it, v) € with H(x) = H+ is given by 
a (small) normal displacement e(x). In |^] we show that 
the corresponding variation 5H(x) = H (x') — H± is given 




FIG. 1: (a) The Schwartz P-surface, a minimal (zero mean 
curvature) surface of simple cubic symmetry. A single unit 
cell is shown, (b) A single unit cell of size L of a CMC surface 
with HL « 2.0, which is close to a simple cubic lattice of 
touching spheres. Figures courtesy of J.T. Hoffman, MSRI 
Scientific Graphics Project, UC Berkeley JhJ. 



to first order in e by 

£ff(x) = - (2Hl - K)e + V a B ■ V s e 
S(x) = -i log (if* 2 - if (x)), 
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where V s and V s 2 are the surface gradient and Lapla- 
cian, respectively, and K(x) is the Gaussian curvature 
at x. Substituting Eq. (Q) in (^), we obtain the linear 
dynamics near a CMC surface with mean curvature 
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Recall the definition of H and K in terms of the local 
principal radii of curvature of the surface R\ and i?2 , 

H = (R^ 1 + Rz 1 )/2, ; K = (R^)" 1 . 

Since 2H 2 - K = (R~ 2 +R~ 2 2 )/2 > 0, uniform precipita- 
tion (e > 0) decreases the curvature, making the surface 
locally more hospitable to deposition, while dissolution 
(e < 0) increases the curvature, favoring further dissolu- 
tion. This observation already suggests that the surface 
is likely to be unstable. 

Because the term in C which couples Gaussian curva- 
ture gradients with the gradient of e is non-Hcrmitian, it 
is more convenient to work with the Hermitian operator 
TL defined via the "gauge transformation" 



so that 
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FIG. 2: (a) The curvature HL versus volume fraction <j> 8 for 
the constant mean curvature extension of the Schwartz P- 
surface with unit cell size L. The minimal (H=0) P-surface 
is represented by the point M, and the endpoints correspond 
to (A) a simple cubic lattice of touching spheres, and (A') 
a simple cubic lattice of touching spherical holes, (b) The 
surface area /L 2 per unit cell of size L versus enclosed solid 
volume fraction <j> 3 for the same simple-cubic CMC family. 
The lower branch is stable under periodic, volume-preserving, 
surface-minimizing dynamics. Figure adapted from Ref. B, 
courtesy of Advances in Chemical Physics. 



The eigenvalues of TL are identical to those of C Jll| , and 
they are bounded from above provided that V(x) < oo. 

For any function 'J'(x) defined on the surface, a lower 
bound on the maximum eigenvalue ujq of TL (consequently 
C) can be established using the inequality 



f T dS^(TL^) 
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where J r dS is the ordinary surface integral. Using the 
trial function ^'b(x) = e s ( x \ we obtain 
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However, since H 2 - K = [(l/i?i) - (l/i? 2 )] 2 / 4 is non- 
negative everywhere on the surface and reaches zero ho- 
mogeneously on r only for a planar surface, we conclude 
that the only marginally stable fixed point of the dynam- 
ics ^ is the plane. There are no other surfaces of any 
kind, porous or non-porous, that are metastable under 
Ostwald ripening. This conclusion applies both to pe- 
riodic and to possible amorphous CMC surfaces. For a 
system with a pore scale of L, we expect this instability 
to manifest itself on a time scale of 



Note that intuitively, the instability is enhanced by an 
increase in the magnitude either of the mean curvature, 
or of the Gaussian curvature (since the latter is on aver- 
age negative for a surface of large genus, such as a typical 
porous surface.) 

For triply periodic CMC surfaces, the periodicity of the 
"potential" V(x) implies that the eigenmodes Vl / k,ra(x) of 
the operator TL (and of C) are surface Bloch functions, 

W* k ,„(x) = Wk, n e 4k ' x £/ k ,„(x) (17) 
t/k,n(x) = (7 k ,„(x + R) (18) 

where R is a lattice vector, k is the crystal momentum 
(confined to the first Brillouin zone), and n is a discrete 
label distinguishing between different branches. Stan- 
dard techniques should allow determination of the spec- 
trum of TL throughout the first Brillouin zone |)| . 

As a concrete example, we have computed the posi- 
tive eigenvalue luq corresponding to the k = eigenmode 
^o(x) = t7o,o(x) for the Schwartz P-surface numerically 
by maximizing the functional on the RHS of Eq. (pl|). 
The details of the maximization procedure will be pre- 
sented elsewhere |)|. We obtain loq « 6.1 for a unit cell of 
size L = 1, whereas the lower bound implied by Eq. ( |l5| ) 
is ojq > 2.9. 

Since the eigenfunction of the largest eigenvalue has 
the same sign (precipitation or dissolution) everywhere 
on the surface, the unstable dynamics of Eq. (^) involves 
transport of solute into (or out of) the pore space from 
the distant reservoir. Therefore, the nature of the dy- 
namics might be different in a closed system, where the 
total amount of the solid component is conserved. 

In a "perfect closed" system, we further regard the 
surface motion as reaction limited, such that the solute 
concentrations near the entire surface are approximately 
equal to the average concentration in the pore c. This 
condition is satisfied at the pore scale L if D ^> KfL/c sa t- 
By integrating Eq. (||) , we see that c must then approach 
a value c oq determined by 
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controlled by the average mean curvature H of the sur- 
face, over a characteristic time scale r pore = Lc sa t/Kf. 



For r > t d 



the solute concentration c is "slaved" to 



j(i/). Conservation of the total solid requires 

4>t = 4>b + c(l - S ) , 



(20) 



where <p s is the solid volume fraction enclosed by the 
surface and <j) t would be the solid volume fraction if none 
of the solid were dissolved. Eliminating c = c cq (_ff ) from 
Eq. ( p0| ) using Eq. (|l^) shows that <p s is very insensitive 
to the value of 
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where we have used Eq. <S). In a "perfect closed" system, 
we take 4> s constant, which is an excellent approximation 
for r > T porc . This dynamics is now given by Eq. (||) 
with the constant H+ replaced by H, whose evolution 
with time is controlled by the constraint of constant <j) s ~ 
{(j>t — c sat )/(l — c sat ). This is surface area minimization 
under the constraint that the volume contained within 
the surface is conserved. 

Consider again periodic surfaces of the Anderson type. 
Figure ||(b) shows the total surface area as a function 
of 4> s along the manifold of CMC surfaces of simple cu- 
bic structure. Although all CMC solutions represent 
fixed points of the dynamics, we anticipate that for a 
given <fi s , only the CMC surface with the lowest area 
can be stable under these dynamics, and only if the sys- 
tem is also required to maintain its spatial periodicity. 
We have confirmed this by direct simulation using Sur- 
face Evolver||, This result is an extension of the 
well-known stability of the minimal surface (point M in 
Fig. ||) to the entire lower branch in Fig. |](b), under peri- 
odic surface area minimization dynamics with conserved 
volume. A corollary of this result is that along this lower 
branch, C has only one unstable k = mode 0, which 
is disallowed by solid conservation under closed system 
dynamics. Thus the stable mode with the longest relax- 
ation time T re iax will control approach to the stable CMC 
surface; dimensional analysis suggests that T rc i ax ~ r ripon . 

Nevertheless, for an extended system this apparent 
stability is compromised by the unstable hydrodynamic 
(k ^ 0) modes in the eigenvalue spectrum Eq. (fL7|), which 
do conserve the overall <p s . It is possible for these un- 
stable modes to be restricted to the topmost (n = 0) 
band and to wavevectors k < tt/£q, such that transport 
of solute over length scales > £q is needed to activate 
these modes. For amorphous surfaces, the band struc- 
ture picture is not appropriate, but we still expect un- 
stable modes to appear above a characteristic length scale 
Iq > L. In either case, it is natural to ask if there is a 
window of time over which the CMC surfaces might be 
observable, before the instability manifests itself. 

During the time r re i a x characterizing the relaxation 
to the CMC surface, the diffusion length over which 
solid transport is possible is £d — V-D^ciax- Using the 
reaction-limited constraint D 3> KfL/c S3j t and T rc i ax ~ 



Tripcn, with T ripen given by Eq. (|16|), we find 
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where the last inequality follows from Eq. (j5|). Since we 
thus expect £d ^> £q, by the time a surface relaxes on 
the pore scale to a CMC surface, the diffusive transport 
between pores will have activated the unstable dynamics, 
taking the system away from this surface. Thus we do not 
expect an intermediate-time window over which CMC 



surfaces would be seen [[13| . 

Using representative values for Kf and er for quartz 
in water |Q, we get for L ~ 100 /zm that r r i pcn ~ 500 
million years, a time scale over which geological systems 
should be regarded as open. However, a similar estimate 
for limestone yields a value of r r i pen ~ 7 years, which is 
quite short on geological time scales, so that these sys- 
tems will correspond more closely to closed systems p5| . 
Of course for real geological systems these processes in- 
volve complicated multicomponent equilibria, and other 
mechanisms such as pressure solution may dominate as 
well [[l6|. Nevertheless, there are probably geological cir- 
cumstances under which surface-tension driven phenom- 
ena will be important. 
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